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Abstract. The theory of incoherent nuclear resonant scattering of synchrotron radiation 
accompanied by absorption or emission of phonons in a crystal lattice is developed. The 
theory is based on the Maxwell's equations and time-dependent quantum mechanics un- 
der the condition of incoherent scattering of radiation by various nuclei. A concept of 
coherence in scattering processes, properties of the synchrotron radiation, and conditions 
of measurement are discussed. We show that employing the monochromator with a narrow 
bandwidth plays a decisive role. The equations for energy dependence and time evolution 
of scattered radiation are derived in detail for the case of the non-split nuclear levels. The 
effect of hyperfine interaction on the time spectra is also considered. We show that the 
Singwi and Sjolander theory of inelastic nuclear absorption, improved by the convolution 
with the instrumental function of the monochromator, is valid for the incoherent nuclear 
scattering in the case of integral over time measurements. 
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1. Introduction 

Nuclear resonant scattering of x rays attracts much interest after the Mossbauer's dis- 
covery of pure elastic nuclear resonant scattering . The recoil energy is transferred 
to an entire crystal lattice, and the frequency of the scattered by nucleus radiation 
remains resonant. This gives an access to studies of coherent phenomena in nuclear 
scattering and stimulates a development of theory of coherent elastic nuclear resonant 
scattering (see the review of theoretical results in Refs. [2,3]. It became clear soon 
after the discovery that, because of the very small energy width of nuclear resonance, 
nuclear resonant scattering gives a good opportunity for studying dynamic of nuclear 
motions in crystals, especially thermal vibration and diffusion [4-6]. However, mea- 
surements of phonon spectra with radioactive sources were extremely difficult due to 
a small cross-section of inelastic processes. 
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Recently the situation has drastically changed owing to the development of the 
synchrotron radiation sources of third generation such as ESRF (France), APS (USA) 
and SPringS (Japan). A high brightness and spectral density of these sources, on the 
one hand, and a development of x-ray monochromator with energy bandwidth less 
than one millielectronvolt, on the other hand, give rise to numerous experimental 
studies of nuclear resonant inelastic absorption. These experiments deal with mea- 
surements of atomic fluorescence, which results from resonant absorption of incident 
photon by nucleus with simultaneous absorption or emission of phonons in solids. 
The review of experimental works and details of experimental methods are discussed 
in Ref. 0. The analysis of the experimental results in the first works is based on 
the theory by Singwi and Sjolander p, which is, strictly speaking, formulated only 
for the particular case of cubic Bravais crystal lattice. The general theory of nuclear 
inelastic absorption for anisotropic crystals is developed in Ref. 

We note that the scattering process is not discussed in these works. Only that 
part of nuclear absorption cross-section, which corresponds to absorption of photon 
by thermally vibrating nucleus is considered. The incident radiation is assumed 
as a monochromatic wave, i.e., the specific properties of the synchrotron radiation, 
the details of measurement and the time dependence of scattered radiation are not 
considered. One can believe that the theory is approximately valid for the nuclear 
resonant inelastic absorption accompanied by the atomic fluorescence, however, being 
not adapted to the real experimental conditions. As for the theory of nuclear resonant 
inelastic scattering (nuclear fluorescence) the problem remained open before recently. 
An attempt of developing the theory of both the nuclear fluorescence and the atomic 
fluorescence, which follow the process of inelastic nuclear excitation, is presented in 
Ref. using methods of quantum electrodynamics and Green functions. However, 
in our opinion, the theory is not complete in the part where the details of experiment 
are discussed. Also the complicated technique does not allow obvious description of 
the scattering processes. The theory is difficult for understanding and requires special 
theoretical background. 

This limitation can be overcome in frames of another approach which, after the 
Laue's work, is widely used for description of various scattering processes. It is 



based on the Maxwell's equations and dielectic properties of medium. In many works 
(see, for example, Ref. p|) nuclear resonant scattering is described by Maxwell's 
equation where a source of scattered radiation is treated as nuclear currents calculated 



by means of the time dependent perturbation theory of quantum mechanics [|rT| . In 
this work we present the theory of incoherent nuclear inelastic resonant scattering 
within this conventional approach. This allows us to take explicitly into account 
properties of the synchrotron radiation and conditions of measurements. 



2. Maxwell's equation and coherence 

For X rays with the energy about 10 keV and more, the electric E(r, t) and magnetic 
B(r, t) fields are strictly connected to each other, and the values of the modulo squared 
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field strength are identical after averaging over a period of oscillations . Therefore it 
is sufficient to consider an equation for any of these fields, or for some vector which 
allows one to determine both of them. In the quantum theory of radiation (see, for 
example, the vector potential A(r, t) in the Coulomb gauge is used. It is chosen 
in such a way that 

V-A = 0, B = VxA, E = — . (1) 

c dt ^ ^ 

The Maxwell's equation for the vector potential can be written as 

"^'Mr^t) - \^A{r,t) = --}{r,t) , (2) 
Ot'^ c 

where j(r,t) is the current acting as a source of electromagnetic radiation, c is the 
speed of light. 



2.1 Coherence as a possibility to observe interference. 

The current j(r, t) results from moving charged particles. If the particles move in- 
dependently of the radiation, they are regarded as free sources of radiation. When 
the motion of particles is forced by incident radiation, we deal with re-radiation or 
scattering. In the latter case we deal with the induced current, which is proportional 
to the radiation field strength. The Maxwell's equation is linear relative to the field 
strength, and the induced currents for the x rays are also linear with high accuracy. 
Therefore the equation obeys the superposition principle for any ensemble of sources 
and scatterers. In other words, any fields and any scattering processes are coherent 
to each other. 

This well known fact of classical optics should be stressed here, because sometimes 
one speaks about a principal incoherence of some sources or scattering processes. In 
reality, the problem of coherence absence is connected with conditions of measure- 
ment. Though the intensity of radiation from various sources or scatterers is always 
positive and must be only added, the interference term in the modulo squared total 
field strength has alternating signs in various time-space regions. As a rule, it strongly 
oscillates in time and space and its average value vanishes. Therefore, if experimen- 
tal conditions allows only a measurement of average value of the interference term, 
the results cannot show coherent phenomena. However, if an experimental technique 
and properties of incident radiation allows one to observe the interference term, the 
processes are coherent. Thus, in order to observe the coherence one needs specific 
experimental conditions with essentially high spatial or temporal resolution of the 
detector. 

For example, consider excitation of nuclei by x rays with extremely broad energy 
range containing all nuclear resonances. In this case, nuclear resonant scattering 
via various nuclear levels creates waves with so large difference in frequencies, that 
it would be impossible to record the interference term even with the best modern 
technique. Therefore, these acts of scattering should be considered as incoherent 
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[P]. However, scattering via various sub-levels of hyperfine structure of a single nu- 
clear level leads to quantum beats in the time spectra with the period of about tens 
nanoseconds. Time resolution of available detectors allows measurements of result- 
ing interference (see the review of experimental results in Ref. [0). Thus, these 
quantum beats is a result of coherent phenomena. 

Similar situations take place in space spectra. For example, if a source is strongly 
localized and a detector has sufficiently high spatial resolution, one can expect to 



record the interference phenomena even with x rays of high energy . In general, 
the concept of coherence can be formulated as a possibility to observe the interference 
term in the intensity of radiation recorded from many sources and/or scatterers. 

There exist particular conditions, when electric fields, resulted from scattering of 
primary radiation by very many nuclei, are added with the identical phase factor. 
Then the total intensity of scattered wave is proportional to A^^ where is a number 
of nuclei. In this case the interference term can significantly exceed the average 
intensity, and the strength of the scattered wave can be comparable with the strength 
of the incident wave. The effective wave field, which excites the nuclei, is not known 
a priory in this case, and it should be calculated self-consistently. These conditions 
are called dynamical or multiple scattering. A conventional approach of dynamical 
theory of scattering supposes a calculation of the field inside the volume of scattering 
medium solving the self-consistent equations. The medium is considered as infinite or 
limited in one direction (lamina- like medium). The fields at both sides of the entrance 
surface of the sample must be identical (boundary conditions). The field at the exit 
surface is considered usually as that measured by a detector, although successive 
propagation of the wave in air may be important in some cases. In the dynamical 
theory, it is sometimes convenient to introduce a continuous spatial distribution of 
the scatterers and the medium polarization similarly to classical optics. A particular 
case of the dynamical theory is the kinematical theory, where the scattered wave 
remains much weaker than the incident wave, but scattering is still spatially coherent 
and depends on the shape of the sample. 

Kinematical scattering is often identified with the first Born approximation al- 
though this is not completely correct. In the first Born approximation the boundary 
problem is not considered at all. The scatterers are assumed to be far away from both 
the source and the detector, whereas the scattering amplitudes are expected to be 
very small. If the experimental conditions allow one to record the interference term, 
the scattering process is coherent. If the interference term can be expected as negli- 
gible, it is sufficient to calculate the intensity of scattered wave from one nucleus. A 
total intensity of radiation becomes proportional to a number of nuclei in a target. In 
this approach, however, the act of scattering should not necessarily be connected to 
a single nucleus. In general, the incident wave illuminates many nuclei. The problem 
of localization or de- localization of excitation is, in fact, identical to the problem of 
phasing or de-phasing of the interference term in the total intensity of the radiation 
scattered by many nuclei. 

In this work, we consider nuclear resonant scattering accompanied by absorption 
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or emission of phonons in a target. The process has relatively small cross-section. 
Therefore scattering events are collected within a large solid angle and are accumu- 
lated after many flashes of incident synchrotron radiation. These conditions do not 
allow us to record the interference between various nuclei. Thus, we can restrict 
ourselves by an analysis of the incoherent scattering. We will use the first Born 
approximation in the problem of scattering by a single nucleus. A summation over 
many nuclei leads to a trivial factor of a number of nuclei in the unit volume of the 
target. The task is solved in three independent steps. Firstly, we declare the proper- 
ties of the incident radiation Aj(r, t) = ejy4j(r,t). Secondly, we calculate the current 
j/j(r,t), of nuclear resonant transitions. Finally, we calculate the scattered radiation 
Af{r,t) = efAf{r,t), at large distance from the nucleus. Here ej ej are the unit 
vectors of polarizations. 



2.2 Temporal and spatial structure of incident wave. 

Properties of synchrotron undulator radiation are well known. They are complicated 
by significant relativistic phenomena. Instead of using them straightforwardly, we 
consider a simplified model, where the current of the source is strictly localized in 
time and space. 

1 1 

V^A{r, t) - ^g^A{r, t) = -- 6{r - ro)<5(t - Q . (3) 

Here Tq is the source position, and tg is the instant of flash. Although a bunch of 
electrons in a storage ring consists of many electrons, the radiation from various 
electrons can be considered as incoherent. Therefore each event of nuclear scattering 
is initiated by incident radiation from a single electron. The width of the synchrotron 
radiation pulse (e.g., about 100 ps at the ESRF, Grenoble) does not enter the problem. 
The solution of eq. is obtained as the Fourier integral 

, , , r du , . r IN exp(iu;|r — rnl/c) 

A(r,i) = /^exp(->.[i-i.l)^^L_L— iO, (4) 

where the particular solution for the monochromatic source is applied as a spherical 
wave. The integral in (^) can be calculated analytically as follows 

«(«-«,-^-r^^ (5) 

47rc|r — ro| 

This result has a simple physical meaning. Instantaneous localized excitation spreads 
in all directions from the point ro and is distributed over the surface of the sphere of 
radius |r — ro|. At the instant t the radius is equal to the distance |r — ro| = c{t — tg) 
which the light covers during the time interval {t — tg). The intensity of the radiation 
is proportional to |r— ro|^^. The total intensity over the sphere area remains constant. 

In further calculations, the integral presentation (^) is more convenient. At large 
distances from the source |ro| ^ |r| we can use the approximation 

^ 47rc|ro| / ^ ^^P(^o^ - iuj[t - - to]) , (6) 
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where to = |ro|/c is the flight time through the distance |ro| from the source to 
the origin of the our coordinate system where the nucleus is located, kg = {uj/c)sq, 
Sq = — ro/|ro|. We use the approximation |r — fqI ~ |ro| + Sor. The expression 
describes the set of plane monochromatic waves with all possible frequencies. 

In reality, an electron in the storage ring is not localized in space but moves with 
the near-light speed. Therefore, the synchrotron radiation is strongly anisotropic. It is 
close to the spherical wave only in the electron coordinate system ||T^, |T6| . At the point 



far from the source and in the laboratory coordinate system, it can be approximated 
by the Gauss beam directed along the speed of motion. Within its central part, 
however, the radiation wave field is close to the spherical wave which approximately 
coincides with the plane wave. Thus, the postulated spatial localization is consistent 
with experimental conditions. In order to estimate the temporal localization of single 
electron radiation flash, we can consider an effective time width as given by frequency 
spectrum of synchrotron radiation under typical experimental conditions and the 
uncertainty principle. The estimated time width turns out to be negligibly small 
compared to the nuclear hfe time. Thus, eq.(P approximately describes the incident 
synchrotron radiation under the typical experimental conditions. 

In experiments on nuclear inelastic scattering the bandwidth of radiation is lim- 
ited by a monochromator. Let the transmission amplitude of the monochromator be 
P{uj — ujq), where c^o is the frequency at the maximum of transmissivity through the 
monochromatoron, while the instrumental function Im{^) = \P{^)f is centered at 
zero argument. The frequency band, provided by the monochromator, is sufficiently 
small to neglect the frequency dependence of the wave vector kg, because the consid- 
ered values of the vector |r| are small and comparable with displacements due to the 
thermal motion of the nucleus. As a result, we arrive to the following approximation 
for the amplitude of the incident radiation 

Ai{v, t) ^ — ■ — ■ / — ex.Y>{-iuj[t -ts- to\)P{uj) . (7) 



47rc|ro| J 2% 

Hereafter kg = (a;o/c)so. 

The frequency integral in eq. (0) describes time evolution of the incident radiation 
flash after passing through the monochromator. It is determined by the response 
function of the monochromator 

P(t)= /^expHc.t)PM. (8) 

According to general physical principles, P{uj) as a function of complex variable u 
can have the poles only in the bottom half-plane. Therefore the response function 
P{t) vanishes for t < 0. In other words, it describes the retarded response. The 
nucleus perceives the radiation flush at the instant t = tg + Iq, which is by to later 
than ts - the instant of the flash at the source . Time to is required for light to travel 
the distance |ro|. The accurate time profile of the flash perceived by the nucleus is 
determined by the monochromator and is described by the function P{t). 
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2.3 Temporal and spatial structure of scattered wave. 

It will be shown below that the current of nuclear transition can be expressed as an 
integral over a limited range of frequencies near the middle frequency cui, 

e/j/i(r,i) ^ j ^ex.p[-i{uJi+uj)t]J{r,u) . (9) 

Therefore we can write the amplitude of the scattered wave at the detector position 
also as a spherical monochromatic waves: 

where ta is the instant of arriving the radiation flash to the detector, and k[ — 

{ui +uj)/c. 

We assume that the distance |ri| between the detector and the scattcrcr is much 
larger than the typical displacement of the nucleus |r| due to thermal vibrations. We 
place again the centre of the scatterer at the origin of the coordinate system and use 
the relation ri — r = |ri|si — r, where Si= ri/|ri|. Then, applying an approximation 
|ri — r| |ri| — s^^r, we arrive to 

Af{ri,U) = J ^^^^^i^^ J c?rexp(-zkir)J(r,u;) , (11) 

Here ki is the scattering vector, ki = fciSi = A;iri/|ri|. Its modulus is ki = Wi/c, 
because we can neglect small frequency u due to the same reason as for the incident 
radiation. The value ti — |ri|/c is the time interval required for light to travel the 
distance |ri| from the nucleus to the detector. 

Thus, the scattered radiation from single nucleus is described by an anisotropic 
spherical wave, which shape is determined by properties of the nucleus. 

3. Quantum current of nuclear transitions. 

The scatterer is the nucleus located in the node of crystal lattice of the target. The 
state of scatterer can be described in terms of the collective state of the centres 
of mass of all nuclei in the crystal \xi) and the internal states of the nucleus 
Consider first the internal nuclear transitions. 

3.1 The Schrodinger equation. 

The current j(r, t) is the quantum mechanical average value of the current density 
operator over the states in presence of the electromagnetic field of radiation. Since 
the notation is used for the stationary states of the nucleus we will denote the 
wave function in the presence of the radiation field by \ipj{t)). 
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We consider a 3-level system with \4>i)fi{t) as the initial state, \(pj)fj(t) as the 
intermediate state and \4>f)ff{t) as the final state of the nucleus. The equation for 
the quasi- stationary states is 

^o|0i) = Ej\(j)j) , fj{t) = exp{-iujjt - 'jjt) , Uj = Ej/h , = Tj/2h . (12) 

Here Hq is the Hamiltonian of the nucleus. We explicitly take into account a decay of 
the excited states with time. For the sake of simplicity, in this section we assume, that 
the hyperfine splitting of the nuclear levels is absent. In the problem of the energy 
spectrum of inelastic nuclear scattering with participation of phonons the hyperfine 
splitting is not essential at all. It infiuences only the temporal characteristics of 
scattering. This question will be analyzed later. 

The ground state of the nucleus \ipiit)) is perturbed by the field of incident radia- 
tion Aj(r, t). The Hamiltonian Hint{t) of interaction of the nucleus with the radiation 
field is small and can be taken into account in frames of perturbation theory. In 
general, the perturbed ground state contains small contributions from other nuclear 
states: 

\m) = mm + (^Am)m + ■■■■ m 

This function is the solution of the Schrodinger equation 

th^^ = [Ho + HUt)Mt) ■ (14) 
A substitution of (|T3p into (0) leads to 

^^-Q^4>i + (^ji^^-Q^4>j + i^-^4>jfj ~ Ho(j)ifi + ajiHo(j)jfj + Hint{t)(f)ifi . (15) 

Since the first and second terms from both sides are equal to each other, we obtain 
f)fi ■ . 

^n^ = {<P,\HUtM^)Mt)fi\t). (16) 

Assuming the adiabatic inclusion of the interaction, i.e. Hint{—oo) = 0, we arrive to 
expression 

= ^ / dt' e:^Y>{iujjit' + '^jt'){(pj\Hint{t')\(pi) ^ (17) 
ih J —oo 

Here We explicitly took into account that 7^ = 0. 



3.2 Induced current as a source of scattered radiation. 

Let us consider the matrix element of the current density operator for the transition of 
the second order from the initial state i into the state / assuming that the transition 
is realized through the intermediate state j which is excited by a radiation transition 
from i into j. The final state / must be considered with zero radiative width, because 
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a finite width describes transition probability and, therefore, has sense only for the 
intermediate state. Applying the Fourier transformation we have 

j^^^"'^^ = / (^^^'^^""^*^^^^^''('^/ Ij'-(k)IV^^W) ' (18) 

where r„()f:) is the coordinate of the vibrated nucleus. Resonant scattering is described 
by the perturbed term in In accord with (|T^ we have 

X Hint{t')\(|)^) exp[i{-Uji + n,){t - t')] . (19) 

The Hamiltonian of interaction is well known. It is convenient to write it in the 
form 



H^nt{t')=-^ J rfr'j„(r'-r„)Ai(r',t' 



=--Jd-'J ^e^'^'(^'--(*'))eJ„(kV.(r,t'), (20) 

where, as above, is the unit vector of polarization of incident radiation. Then, the 
projection of the current density on the unit vector of polarization of the scattered 
wave may be written in the form 

e;j^,(r,t)= f dt' j dv' Mf,{v,t-v\t')A{v\t'), (21) 

J —CO J 

where we introduce the scattering matrix 



X exp[^(-c^}"^ + i^j){t - t')]e-^kr„(t)g-ik'r„(i') _ ^22) 

Hereafter the frequencies cjj-"'' related to the nuclear transitions are marked by the 
index (n) and we use notation 

iV;,(k,k') = (0;|e;L(k)|0,)(0,|eJ„(k')|0.). (23) 



Expression ( P^ ) for the scattering matrix contains the coordinates of the nucleus at 
the instant of excitation t' and at the instant of relaxation t. 

The nuclear coordinates depend on the time owing to thermal vibrations. As 
known, the thermal vibrations of atoms in a crystal lattice exhibit distinct quan- 
tum behaviour, thus they must be described as absorption and emission of phonons. 
Therefore we will consider the coordinate r„ as the operator and introduce the wave 
function of phonons. 

m 

=exp(ia;5?^t)L/i(k, k', t~f). (24) 
9 



Here index (p) shows that the frequency is related to the phonon system. The sum 
over intermediate states of the phonon system is required in order to describe an 
evolution of the system from the instant t' into the instant t. The function Lfi 
depends only on the time interval and is defined by expression 

Lf,{k, k', t) = J2{xf \e-'''''"\Xm){Xm le-^'^'^lx.) exp(-2c.St) . (25) 

m 

As a result, the scattering matrix takes the form 
Mfiir, t- r', exp[2(-cu, + «7)(t - t')] 

X / ^e^'^^^''^'A^/.(k, k')L,,(k, k', t-t'), (26) 



(2vr)^ 

where we introduce the total transition frequency oofi = uj^^i + cuj^^ the resonance 

frequency ujr = ooj^^ and omit the index of the width of the excited state 7 = 7j. 
Expression (|26| ) is close to eq.(25) of Ref. 0. However, here it is derived for 3- 
level system within standard quantum mechanics approach. A summation over the 
intermediate states of the nucleus, as made in Ref. 0, does not have sense because 
the resonance can be realized only for one transition. 

Substituting eqs. (|26| ) and (|^ into eq. ( pi]) and making identical transformations 
we obtain the explicit form for the function J(r, u) in eq. @ as follows 

J(r, uj)=PH ^-^^ J ^ exp(.kr) Ar,.(k, -ko) 

xLfi{k,-ko,uJo + uj-uJr + i'y) (27) 

Here loi = uq —ujfi is the middle frequency of the scattered wave. We also introduced 
the function 

Lfi{k,k',uj) = dtexp{iujt)Lfi{k,k',t) . (28) 



4. Properties of scattered radiation. 



Substituting (|27| ) in ( |TI1) we obtain the vector potential of the scattered wave as 
follows 

(47r)^c'^ft|ro||ri| 

X [ ^exp{-iujt)Lfi{ki,-ko,uJo + 00 - uJr + il)P{(^) ■ (29) 
J In 

Here t = t^t — ~ — ti is the time delay of radiation in the monochromator and 
in the nucleus, because is the time instant of arrival of the scattered flash to 
the detector, whereas {tg + to + ti) have sense of the time instant of arrival of the 



flash propagating in vacuum along the identical trajectory. The definition (P8|) states 
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that the function L/j(ki, — ko,ix;) has the poles only in the bottom half-plane of the 
complex variable uj. Therefore the strength of the scattered radiation as function 
of the delay time t vanishes for t < 0. As a result, we obtain general expression 
for intensity of scattered radiation as a function of the delay time t and the mean 
frequency ujq of the monochromator transmissivity as 

2 



I{uJo,t) = IqIn 



duo 
2^ 



exp(— ic<jt)I//j(ki, — ko, ojq + u — uJr + i'^)P{(jj] 



where 



Ary,(ki,-ko 



(30) 



(31) 



' (47r)4c8^Vgr? 

In measurements of nuclear inelastic scattering, the signal is accumulated over many 
pulses of synchrotron radiation and resulted from scattering by many nuclei. There- 
fore the above expression must be summed over the final states (index /) and averaged 
over the initial states (index i) of the phonon system. We will denote such intensity 
as I{uo,t) and call the average intensity. 



4.1 Energy dependence of scattered radiation. 

It is instructive to analyze first the complete energy spectrum of the scattered radi- 
ation when the measurement is performed during the total delay time interval from 
zero to infinity and the integral intensity is 



I{uJo, t)dt , 



(32) 



Since I{t) = at t < 0, the integration limits can be formally replaced by (— cxd, oo). 
Taking into account that expression ( |5D| ) is the modulo squared of the Fourier integral 
over frequency, we use the Parseval's theorem and obtain straightforwardly 



hntiujo) = IqIn I ^ Lfi(ki, -ko, uq - Ur + uj + ii) hiiuj 



(33) 



Thus, the time-integrated intensity of scattered radiation as a function of the middle 
frequency ujq of incident radiation is a convolution of the energy spectrum of inelastic 
scattering for the case of monochromatic incident wave with the instrumental function 
as the energy spectrum of the monochromator Im{^^) = |-P('^)|^- 

In order to calculate the averaged intensity Iint{uJo) we have to average the function 

Lfi(ki, — ko,u; + ^7) . Substituting ([25| ) in (pSf), we obtain 



L/,(ki,-ko,u; + i7) = -J2 



{Xf |e" 



-jkir„ 



Xm) {Xr. 



Xi) 



[UJ — UJ. 



ip) 



(34) 



The square modulus of the right-hand side must be summarized over the final states 
that leads to the expression 



E 

f,m,m' 



(X/|e 



\Xm){X r, 



o«knr„ 



\Xi) {XiW 



-«kor„ 



\Xm' ){x r. 



e 



«kir„ 



\Xf) 



(^ - + il) 



a; — u. 



(p) 



(35) 
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Here the sum over / can be calculated independently 



{Xm' W^'"'^Xf){Xf |e '^'""^Wm) = {Xm'WXm) = ^mm' • (36) 



Then the dependence on the direction of the scattered wave vanishes and, instead of 
the triple sum, we obtain only a sum over m . 

Replacing index m ^ f, we may present the averaged intensity in the form 

hnti^^o) = IqIn J ^Issi^o -ujr + uj) /m(c^) , (37) 
where 

T ( \ /X-\f n 1 ^- ^^ ^ 9i\{Xf\ exp(zkor„)|xi)P . . 

IssM = {2^\Lf^{kl,-]^o,uJ + ^-f) }t = 2^ -73 ; • (38) 

/ ij (^//-w)^ + 7^ 

Here (• • ■)t means the thermal averaging over the initial states, and Qi is the weight 
of z-th initial state of the phonon system. 

Equation ( |57| ) is a standard recipe of taking into account the instrumental func- 
tion (see, for example, Refs. p, |18[). Eq. ( pS]) for the integral intensity of scattered 
radiation coincides with the known expression for cross-section of nuclear resonant 
absorption, used by Singwi and Sjolander (see eq.(l) in Ref. 0). An important 
property of this formula is that absorption or emission of phonons takes place simul- 
taneously with resonant absorption of incident photon by nucleus, i.e., at the instant 
of nuclear excitation. The intensity of scattered radiation depends only on the wave 
vector ko, whereas the wave vector ki does not influence the energy spectrum. As 
follows from eq . (^61) , this feature is a logical sequence of the procedure of summation 
over all the possible final states of the phonon system. This treatment corresponds 
to measurement of scattered quanta independently of their energy and directions of 
propagation. If an experimental setup would contain additional elements allowing se- 
lection of scattered radiation with a particular energy, or if scattered radiation would 
be accepted within essentially small solid angle, we have to use instead of eq.(pU[) 
some more complicated expressions with an explicit integration over frequency lvi. 
Besides that, we would have to use the theory of coherent scattering, i.e., to take 
into account interference between various nuclei. In general, energy dependence of 
radiation scattered into some particular directions may then differ from the averaged 
value. 

Thus, under the conditions of entire averaging over angle and energy of scattered 
radiation, the energy dependence of scattered radiation coincide entirely with the 
energy dependence of nuclear inelastic absorption. This conclusion agrees with the 



experimental results of Ref. ^7j. One of the reason of this coincidence is that 
spontaneous emission of photon during de-excitation of a single nucleus does not 
depend on the resonance condition of preceding nuclear excitation. In particular, it 
does not depend on the phonon state, which has determined conditions of nuclear 
excitation. That is why, the energy spectrum of incoherent nuclear inelastic scattering 
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does not differ from tliose of nuclear inelastic absorption. At the first step, a nucleus is 
excited by an incident quantum with some energy transfer to absorption or emission 
of phonons. The following evolution of excited nucleus may vary, but it does not 
influence the dependence of recorded radiation on the energy of incident radiation. 

Following Ref. 0, it is convenient to represent the function Iss{^) as integral 
over time. Below we use slightly different, than in Ref. |^, way of derivation which 
is more useful for us. Let us write the expression (|38|) as follows 



Issiuj) = Y.9^\{Xf\ exp(^kor„)|x.)r^(a;5-? - u) , (39) 
where 

R{u;) = , / = / dtexp{tut) R{t) . (40) 
Substituting the Fourier integral (|40| ) in eq.(|39D we obtain 

IssH = J dtexp{-tujt) F(ko,t) /?(t) , (41) 

where the time correlation function F{hQ,t) depends on the lattice vibrations and is 
determined by 

F{^o,t)=J29iY.(Xi\ exp(-«kor„)|x/)(x/l exp(ikor„)|xi) exp(ia;j?^t) 

i f 

=Yl9i{Xi\ exp(-zkor„(0)) exp(zkor„(t))|xi) 

i 

= (exp(-ikor„(0)) exp(ikor„(t)))T. (42) 
The function R{t) may be calculated in explicit form 

R(t) = ?^^. (43) 

As a result, we arrive to expression 

lintM = IqIn J dtexp{-t[uJo - cj,]t)F(ko, t) R{t) /A/(t) , (44) 
where 

hiit) = f ^exp{-iLut)lM{uj). (45) 



In Refs. p, similar expression was used in order to obtain phonon density 
of states from energy dependence of nuclear inelastic absorption. Our results show 
that, if incoherent nuclear inelastic scattering is ideally integrated over time, its en- 
ergy dependence is identical to that of nuclear inelastic absorption. Strictly speaking, 
an integration over infinite time interval cannot be reached in real experiments, be- 
cause during a short period after synchrotron radiation flash an acquisition system 
is gated against strong electronic scattering. In order to analyze possible influence of 
this gating on energy spectra of inelastic scattering, we discuss time dependence of 
scattered radiation intensity. 
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4.2 Time dependence of scattered radiation. 

Eq. (^) for intensity of scattered radiation can be presented as follows 

I{uJo,t) = IoIn \Qfi{ooo,t)\\ (46) 
where 

Qfi{uJo,t) = J ^ exp{-iujt)Lfi{ki, -^0,1^0 - Ur + uj + i-f)P{uj) . (47) 
Substituting ( p8|) we write this expression in the form 

Qfi{uo,t) = f dt'exp{i[uo-uJr + h]t')Lfi(ki,-ko,t')P{t-t'). (48) 





Here we use eq.(|D and that P(t) = for t < 0. As follows from eq. (|48|), the intensity 
vanishes at t = 0. Making a replacement of the integration variable t' —>■ t — t' and 
using eq. (p5|) we obtain 

m 

X [' dt' exp{i[iUr + J^l - uJo]t') e^'' P{t') . (49) 



General evaluation of this expression requires an explicit form of the response function 
of the monochromator. We will consider the case where the energy width of the 
instrumental function is much larger than that of nuclear resonance. Therefore, the 
response function P{t) decreases significantly on the characteristic time tp, which is 
much less than the life time of excited nuclear state to = 2/7. Then the function 
exp(7t') under the integral can be approximated by unity. It is evident that within 
the time interval < t < tp the integral in eq. (^) increases with time and it 
becomes constant when t ^ tp. In this latter region we can replace the upper limit 
by 00. Omitting also a constant phase factor which does not infiuence the intensity, 
we obtain 

QfiiiUo,t) = e-^'Y.iXf \e-'''''"\Xm){Xm \e^''''-\x^)e~'^-'PiiOr + - c^o) . (50) 

m 

Repeating derivation as discussed above, we arrive to the following expression for the 
intensity summarized over the final states and averaged over the initial states 



I{uJo,t) = Ioexp{-t/to)lNj29i\{Xf\ exp(ikor„)|xi)|^/M(t^r + - i^o) ■ (51) 

We obtain that within the time region t ^ tp the intensity decreases in an expected 
way, i.e. exponentially with the life time of the excited nuclear state. 

Eq. (^) is similar to eq. (0), derived above for the case of nuclear inelastic 
absorption. Therefore, following the same approach, we obtain 

I{ooo,t) = Jo exp{-t/to)lN J dt' exp(-z[cuo - u;.]t')^(ko, t')4/(i') • (52) 
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If the life time of excited nuclear state is, as we assume here, much longer that 
the duration of the flash of incident radiation provided by the monochromator, we 
can neglect time dependence of the function R{t) in eq. and replace R{t) by 

R{0). Then the integrated over time intensity (^) coincides with the approximate 
expression (Q), obtained by other way. 

Time dependence of R(t) in eq.(0) allows one to study correct behaviour of the 
delayed radiation for very small delay time inside the flash provided by monochroma- 
tor. Within this time interval the intensity of scattered radiation increases from zero 
to a certain maximal value. When the flash is over, the major effect is the delayed 
re-emission of excited nuclear state which decays exponentially with the time. 



5. Effect of hyperfine structure of nuclear levels. 

In this section we consider more complicated case, where nuclear levels of both ground 
and excited states are split by hyperfine interaction. In eqs. (JT^), ( P^D and further 
down we have to replace uoi by Ui + Au:^\ Uj by Uj + Aujp \ uj by uuj + Aw^^ where 
indexes a, /3 and 7 denote various sub levels of hyperfine structure for the initial, 
intermediate, and final states, respectively. Accordingly, nuclear wave functions have 
indexes \4>ia) and so on. Eq. p^) for the scattering matrix becomes 



M,.(r,t;r^t> ^— -^-'J^'-- j Nj^^k' ,t - t') 

X exp[z(-Cu]") + t^j){t - t')]e-^kr„Wg-ik'r„(t') ^ (53) 



where co'/7,ia = ujf + /^ujij^ ~uji~ ^uj^\ and the nuclear function A^Jf (k, k', t) acquire 
the time dependence 

A^J°(k,k',)f:) =^(0/-^|ejj„(k)|0j^)(0j^|eij„(k')|0ia)exp(iAu;j^,i<:,t) , (54) 

where Aujf^^ia = Aujj^'' — Ac<j^*\ 

Repeating the derivation discussed above, it is easy to obtain the general expres- 
sion for the time dependence of intensity as 

I{uJo,t) = Iq^q^ |Q]"(a;o,t) , (55) 
where is the statistical weight of the state (m), and 

Q}r(^o, t)= 1^ dt' exp(^[a;o - + ^l]t') A^/f (ki, -ko, t') 

xLf,{k,,-ko,t')P{t-t'). (56) 

We note that hyperfine splitting can significantly exceed the width of the nuclear 
resonance. Nevertheless, the characteristic time ~ l/Ao; of nuclear function 
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variation is much larger than the characteristic time tp of the response function of 
the monochromator. Therefore at t ^ tp we can replace the lower limit of integration 
by —oo and the function, exp(— 7^') Nj^(ki, —kg, t') by its value at the upper limit t. 
As a result, instead of (|52D , we obtain the generalized equation 



J(cuo,t) = Ioexp{-t/to)lN{t) J dt' exp{-t[uJo - uJr]t')F{ko,t')iM{t') , (57) 
where 

/iv(t)=E^"|iV;r(ki,-ko,t)'. (58) 
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To obtain the intensity /Ar(t) integrated over the angle of exit of scattered photons, it 
is necessary to integrate over directions of the vector ki. The function lN{t) does not 
depend on the resonance condition, in particular, on the value {ujq — Ur). Therefore 
it does not influence the energy spectrum of incoherent nuclear inelastic scattering. 



6. Conclusion 



We present the theoretical analysis of energy and time dependence of incoherent 
nuclear inelastic resonant scattering of synchrotron radiation (SR) accompanied by 
absorption and emission of phonons in crystal lattice. The theory is based on the 
Maxwell's equations and the time dependent perturbation theory of the quantum 
mechanics. This allows us to analyze explicitly physical nature of the considered 
processes. 

In order to treat the case of incoherent scattering, we assume that, though the 
scattering process is coherent in general, the conditions of measurements do not allow 
experimental observation of the interference term in the intensity of scattered radia- 
tion due to its smearing out during averaging over various parameters. This concept 
of incoherent scattering significantly simplify calculations, because it avoids calcula- 
tions of the interference term which would in any event vanish during averaging. Our 
results show that if measurements are performed without energy or angular analysis 
of scattered radiation, the dependence of its intensity on the energy of incident radi- 
ation is identical to that of nuclear inelastic absorption. Formally, this is caused by 
averaging of scattering intensity over final states of phonon system. 

The measured spectra are influenced only by processes of absorption and emission 
of phonons only during the nuclear excitation by the SR flash. We assume that x 
rays emitted by various electrons in a storage ring are incoherent, and radiation from 
single electron can be considered as instantaneous. This assumption is based on the 
experimental conditions, where the intensity of scattered radiation is accumulated 
after many SR flashes, whereas only one photon is recorded after each SR flash. We 
show that the possibility to measure the density of phonon states under these condi- 
tions arises from the property of monochromator to delay the initially instantaneous 
SR flash up to the time interval which significantly exceeds the period of nuclear 
vibrations. It is sufficient to describe the case of incoherent scattering within the 
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first Born approximation. Then the calculation is reduced to thermal averaging of 
the phase factor depending on the nuclear coordinates at the instant of photon both 
absorption and emission. Within this approach, the time correlation function can be 
calculated with a proper account for multi-phonon processes, whereas treating the 
multi-phonon processes in the case of coherent scattering is significantly more com- 
plicated problem. 

* * * 

The work is supported by the Russian Foundation for the Basic Research (project 
N. 01-02-16508). 
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